Predicting enzyme class from protein structure using Bayesian classification

Predicting enzyme class from protein structure parameters is a challenging problem in protein analysis. We developed a method to predict enzyme class that combines the strengths of statistical and data-mining methods. This method has a strong mathematical foundation and is simple to implement, achieving an accuracy of 45%. A comparison with the methods found in the literature designed to predict enzyme class showed that our method outperforms the existing methods.

Building multiple sequence alignments with a flavor of HSSP alignments

Homology-derived secondary structure of proteins (HSSP) is a well-known database of multiple sequence alignments (MSAs) which merges information of protein sequences and their three-dimensional structures. It is available for all proteins whose structure is deposited in the PDB. It is also used by STING and (Java)Protein Dossier to calculate and present relative entropy as a measure of the degree of conservation for each residue of proteins whose structure has been solved and deposited in the PDB. However, if the STING and (Java)Protein Dossier are to provide support for analysis of protein structures modeled in computers or being experimentally solved but not yet deposited in the PDB, then we need a new method for building alignments having a flavor of HSSP alignments (myMSAr). The present study describes a new method and its corresponding databank (SH2QS--database of sequences homologue to the query [structure-having] sequence). Our main interest in making myMSAr was to measure the degree of residue conservation for a given query sequence, regardless of whether it has a corresponding structure deposited in the PDB. In this study, we compare the measurement of residue conservation provided by corresponding alignments produced by HSSP and SH2QS. As a case study, we also present two biologically relevant examples, the first one highlighting the equivalence of analysis of the degree of residue conservation by using HSSP or SH2QS alignments, and the second one presenting the degree of residue conservation for a structure modeled in a computer, which , as a consequence, does not have an alignment reported by HSSP.

PDB-Metrics: a web tool for exploring the PDB contents

PDB-Metrics (http://sms.cbi.cnptia.embrapa.br/SMS/pdb_metrics/index.html) is a component of the Diamond STING suite of programs for the analysis of protein sequence, structure and function. It summarizes the characteristics of the collection of protein structure descriptions deposited in the Protein Data Bank (PDB) and provides a Web interface to search and browse the PDB, using a variety of alternative criteria. PDB-Metrics is a powerful tool for bioinformaticians to examine the data span in the PDB from several perspectives. Although other Web sites offer some similar resources to explore the PDB contents, PDB-Metrics is among those with the most complete set of such facilities, integrated into a single Web site. This program has been developed using SQLite, a C library that provides all the query facilities of a database management system

STING Contacts: a web-based application for identification and analysis of amino acid contacts within protein structure and across protein interfaces.

UNLABELLED: Amino acid contacts in terms of atomic interactions are essential factors to be considered in the analysis of the structure of a protein and its complexes. Consequently, molecular biologists do require specific tools for the identification and visualization of all such contacts. Graphical contacts (GC) and interface forming residue graphical contacts (IFRgc) presented here, calculate atomic contacts among amino acids based on a table of predefined pairs of the atom types and their distances, and then display them using number of different forms. The inventory of currently listed contact types by GC and IFRgc include hydrogen bonds (in nine different flavors), hydrophobic interactions, charge-charge interactions, aromatic stacking and disulfide bonds. Such extensive catalog of the interactions, representing the forces that govern protein folding, stability and binding, is the key feature of these two applications. GC and IFRgc are part of STING Millennium Suite. AVAILABILITY: http://sms.cbi.cnptia.embrapa.br/SMS, http://trantor.bioc.columbia.edu/SMS, http://mirrors.rcsb.org//SMS, http://www.es.embnet.org/SMS and http://www.ar.embnet.org/SMS (Options: Graphical Contacts and IFR Graphical Contacts).

Defining 3D residue environment in protein structures using SCORPION and FORMIGA.

SUMMARY: Two web-based applications to analyze amino acids three-dimensional (3D) local environment within protein structures-SCORPION and FORMIGA-are presented. SCORPION and FORMIGA produce a graphical presentation for simple statistical data showing the frequency of residue occurrence within a given sphere (defined here as the 3D contacts). The center of that sphere is placed at the Calpha and at the last heavy atom in the side chain of the selected amino acid. Further depth of detail is given in terms of a secondary structure to which the profiled amino acid belongs. Results obtained with those two applications are relevant for estimating the importance of the amino acid 3D local environment for protein folding and stability. Effectively, SCORPION and FORMIGA construct knowledge-based force fields. The difference between SCORPION and FORMIGA is in that the latter operates on protein interfaces, while the former only functions for a single protein chain. Both applications are implemented as stand-alone components of STING Millennium Suite. AVAILABILITY: http://sms.cbi.cnptia.embrapa.br/SMS, http://trantor.bioc.columbia.edu/SMS, http://mirrors.rcsb.org/SMS, http://www.es.embnet.org/SMS and http://www.ar.embnet.org/SMS. [options: Scorpion, Formiga]

ConSSeq: a web-based application for analysis of amino acid conservation based on HSSP database and within context of structure.

SUMMARY: A web-based application to analyze protein amino acids conservation-Consensus Sequence (ConSSeq) is presented. ConSSeq graphically represents information about amino acid conservation based on sequence alignments reported in homology-derived structures of proteins. Beyond the relative entropy for each position in the alignment, ConSSeq also presents the consensus sequence and information about the amino acids, which are predominant at each position of the alignment. ConSSeq is part of the STING Millennium Suite and is implemented as a Java Applet. AVAILABILITY: http://sms.cbi.cnptia.embrapa.br/SMS/STINGm/consseq/, http://trantor.bioc.columbia.edu/SMS/STINGm/consseq/, http://mirrors.rcsb.org//SMS/STINGm/consseq/, http://www.es.embnet.org/SMS/STINGm/consseq/ and http://www.ar.embnet.org/SMS/STINGm/consseq/

The X-ray structure of a recombinant major urinary protein at 1.75 A resolution. A comparative study of X-ray and NMR-derived structures.

Major urinary proteins belong to the lipocalin family and are present in the urine of rodents as an ensemble of isoforms with pheromonal activity. The crystal structure of a recombinant mouse MUP (rMUP) was solved by the molecular-replacement technique and refined to an R factor and R(free) of 20 and 26.5%, respectively, at 1.75 A resolution. The structure was compared with an NMR model and with a crystallographic structure of the wild-type form of the protein. The crystal structures determined in different space groups present significantly smaller conformational differences amongst themselves than in comparison with NMR models. Some, but not all, of the conformational differences between the crystal and solution structures can be explained by the influence of crystallographic contacts. Most of the differences between the NMR and X-ray structures were found in the N-terminus and loop regions. A number of side chains lining the hydrophobic pocket of the molecule are more tightly packed in the NMR structure than in the crystallographic model. Surprisingly, clear and continuous electron density for a ligand was observed inside the hydrophobic pocket of this recombinant protein. Conformation of the ligand modelled inside the density is coherent with the results of recent NMR experiments.

The high resolution crystal structure of yeast hexokinase PII with the correct primary sequence provides new insights into its mechanism of action.

Hexokinase is the first enzyme in the glycolytic pathway, catalyzing the transfer of a phosphoryl group from ATP to glucose to form glucose 6-phosphate and ADP. Two yeast hexokinase isozymes are known, namely PI and PII. The crystal structure of yeast hexokinase PII from Saccharomyces cerevisiae without substrate or competitive inhibitor is determined and refined in a tetragonal crystal form at 2.2-A resolution. The folding of the peptide chain is very similar to that of Schistosoma mansoni and previous yeast hexokinase models despite only 30% sequence identity between them. Distinct differences in conformation are found that account for the absence of glucose in the binding site. Comparison of the current model with S. mansoni and yeast hexokinase PI structures both complexed with glucose shows in atomic detail the rigid body domain closure and specific loop movements as glucose binds. A hydrophobic channel formed by strictly conserved hydrophobic residues in the small domain of the hexokinase is identified. The channel's mouth is close to the active site and passes through the small domain to its surface. The possible role of the observed channel in proton transfer is discussed.

Crystallization and preliminary crystal analysis of yeast hexokinase PI and PII.

Hexokinase is the prime enzyme of the Embden-Meyerhof pathway and is responsible for the first stage of energy conversion. It catalyzes the transfer of a phosphate to glucose to form glucose-6-phosphate. Yeast hexokinase PII is also known to play an important role in glucose signal transduction. Crystals of yeast hexokinase isoforms PI and PII were obtained by vapour-diffusion techniques using the hanging-drop method. Isoform PI crystals belong to the space group P2(1)2(1)2(1), with unit-cell parameters a = 62.12, b = 78.87, c = 144.74 A. Unit-cell parameters for isoform PII crystals are a = b = 142.81, c = 58.46 A and the space group is I4. Synchrotron diffraction data have been collected to 2.2 A resolution from the isoform PII crystal, whereas isoform PI diffracted to 3.1 A.

Crystallization and preliminary X-ray diffraction studies of piratoxin II, a phospholipase A2 isolated from the venom of Bothrops pirajai.

The phospholipases A2 (PLA2, E.C. 3.1.1.4, phosphatide sn2 acylhydrolases) are the major components of the venom of several snakes. They are responsible for several important pharmacological effects observed in ophidian incidents. PLA2 piratoxin II from Bothrops pirajai has been crystallized by the vapour-diffusion technique. X-ray diffraction data have been collected to 2.04 A resolution (90.2% complete, Rmerge = 0.070). The space group is P21 and the cell parameters are a = 46.19, b = 60.36, c = 58.74 A and beta = 96.05 degrees. The structure has been solved by molecular replacement using the crystallographic structure of PLA2 from Bothrops asper (PDB code 1CLP) as a search model.